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We introduce a measure of both quantum as well as classical correlations in a quantum state, the 
entanglement of purification. We show that the (regularized) entanglement of purification is equal 
to the entanglement cost of creating a state p asymptotically from maximally entangled states, with 
negligible communication. We prove that the classical mutual information and the quantum mutual 
information divided by two are lower bounds for the regularized entanglement of purification. We 
present numerical results of the entanglement of purification for Werner states in 7^2 ® Hi ■ 



I. INTRODUCTION 



The theory of quantum entanglement aims at quantifying and characterizing uniquely quantum correlations. It 
does so by analyzing how entangled quantum states can be processed and transformed by quantum operations. A 
crucial role in the theory is played by the class of Local Operations and Classical Communication (LOCC), since 
. quantum entanglement is non-increasing under these operations. Indeed, by considering this class of operations we 
^ ■ are able to neatly distinguish between the quantum entanglement and the classical correlations that are present in 
.i-j. . the quantum state. 

■ Given the success of this theory, we may be daring enough to ask whether we can similarly construct a theory 
Cn| ' of purely classical correlations in quantum states and their behavior under local or nonlocal processing. At first 

' sight, such an effort seems doomed to fail since merely local actions can convert quantum entanglement into classical 
correlations. Namely, Alice and Bob who possess an entangled state = J2i VKWi) ® Ih) with Schmidt coefficients 
, Xi can, by local measurements, obtain a joint probability distribution with mutual information equal to H{X). Thus 

■ it does not seem possible to separate the classical correlations from the entanglement if we try to do this in an 
Oh, operational way. Note that it may be possible to separate quantum and classical correlations in a nonoperational way, 

• see for example Ref. or The drawback of such an approach is that no connection is made to the dynamical 
Ch I processing of quantum information, which is precisely what has made the theory of quantum entanglement so elegant 
and innovative. An operational approach to the quantification of quantum and classical correlations was recently 
formulated in Ref. 

In this paper we propose to treat quantum entanglement and classical correlation in a unified framework, namely we 
express both correlations in units of entanglement. Such a theory of 'all' correlations may have potential applications 
outside quantum information theory as well. Researchers have started to look at entanglement properties of many- 
^ particle systems for example at (quantum) phase transitions (see for example Ref. [Q and references therein) . Instead 
^ of considering the entanglement of formation in these studies, one may choose to look at the behavior of a complete 
correlation measure. In this paper we introduce such a measure, called the entanglement of purification. We would like 
to emphasize that our correlation measure is not an entanglement measure, but a measure of correlations expressed 
in terms of the entanglement of a pure state. 

It has been the experience in (quantum) information theory that questions in the asymptotic approximate regime 
are easier to answer than exact non-asymptotic queries. Thus we ask how to create a bipartite quantum state p in 
the asymptotic regime, allowing approximation, from an initial supply of EPR-pairs by means of local operations and 
asymptotically vanishing communication. This latter class of operations will be denoted as LOq (Local Operations 
with o(n) communication in the asymptotic regime) versus the class LO for strictly Local Operations. We can properly 
define this formation cost i^LOq as follows: 

{777 
-\3Clo,, i?(/:Lo,(|*-)(*-P"),P^") <e}. (1) 

Here l^*-) is the singlet state in 7^2 (STi.2 and Cloq is a local superoperator using o{n) quantum communication. D 
is the Bures distance D{p,p') — 2-^/1 — F{p, p') and the square-root-fidelity is defined as F{p,p') = Tr(^pi/2p'pi/2) 
. We could have allowed classical instead of quantum communication in our definition, -our results will not depend 
on this choice-, so we may as well call all communication quantum communication. 
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Before we consider this entanglement cost for mixed states, we observe that by aUowing asymptotically vanishing 
communication, we have preserved the interconvertibility result for pure states ||6|. This is due to the fact that both 
the process of entanglement dilution as well as entanglement concentration can be accomplished with no more than 
asymptotically vanishing amount of communication, see Ref. [Q. 

We see that the cost ELOq{p) of creating the state p is defined analogously to the entanglement cost Ec{p) [D, 
1^, with the restriction that Alice and Bob can only do a negligible amount of communication. It is immediate that 
EhOqip) ^ill ill general be larger than Ec{p). In particular, for a separable density matrix Ec{p) = whereas we will 
show that for any correlated (i.e. not of the form pab = Pa® Pb) density matrix E^oqip) > 0. The entanglement 
cost Ec was found |^ to be equal to 

EM - l.m (2) 

n — >oo fi 

where Ef{p) is the entanglement of formation We will similarly find an expression for E^og 

ELOq = lim ' ^ E^{p), (3) 

where Ep{p) is a new quantity, the entanglement of purification of p. 

Our paper is organized in the following manner. We start by defining the entanglement of purification and deriving 
some basic properties of this new function, such as continuity and monotonicity under local operations. We will 
relate the entanglement of purification to the problem of minimizing the entropy of a state under a local TCP (Trace- 
preserving Completely Positive) map. With these tools in hand, we can prove our main result, Theorem ^ Then we 
spend some time proving the mutual information lower bounds for E^oqip)- We also compare our correlation measure 
with the induced Holevo correlation measures Ca/b that were introduced in Ref. [Q. We prove that for Bell-diagonal 
states the correlation measure Ga is equal to the classical capacity of the related 1-qubit Pauli channel. At the end 
of the paper we present our numerical results for Ep{p) where p is a Werner state on 7^2 ®> H2. The proofs of the 
lemmas and theorems in this paper are all fairly straightforward and use many basic properties of entropy and mutual 
information (concavity, subadditivity of entropy, nonincrease of mutual information under local actions etc.). 

II. ENTANGLEMENT OF PURIFICATION 

We define the entanglement of purification: 

Definition 1 Let p be a bipartite density matrix on T-La ®'Hb- Let \ip) G H-aA' ® Ti-BB'- The entanglement of 
purification Ep(p) is defined as 

Ep{p)^ min £;/(|V')(^|), (4) 

V':Tr^'B'IV'>(V'l=P 



where Ef{\il)){il)\) is the entanglement of \'ip) which is equal to the von Neumann entropy S{(Jbb') = — Trcrss' \ogaBB' 
where obb' = '^t:aa' Let {Ai, be the eigenvalues and eigenvectors of pab- The "standard purification" of 

p is defined as 



=Y.V\H^)AB®^A'\i)B'. (5) 
i 

Every purification of p can be written as I?/)) — {Tab ® Ua'B')\'4's) for some unitary operator Ua'B' on A' and B' . 
Therefore, Eq. (^) can be rephrased as: 

Ep{p) = min EHIab ® Ua'B')\A){AKIab ® Ua'B')^)- (6) 

= min S {Tl AA' {I AB ®U A' B')Hs){i^s\{lAB <E)U A' B')^)) 
Ua'B' 

= min 5 ((/b«)AbO(mbb' (/?))), (7) 
where we have taken the trace over A and A' to obtain Eq. (Q), 
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l^BB'ip) = Tr^A'|'0s)(V's| , 



(8) 



and Kb'{v) = Tt:a'Ua'B'{vb' ® |0)(0|^')J7^,^,. The minimization in Eq. is over all possible TCP maps Kb' since 
every TCP map can be implemented by performing a unitary transformation on the system and some ancilla and 
tracing over the ancilla. Note that the minimizations over Ua'B' and Ab' are equivalent. Equations and provide 
two different formulations of the same minimization. Conceptually the first formulation is based on purifications of p 
and variation over Ua'B'- The second formulation is based on extensions of p, (JABB'^ such that Ttb'Ctabb' = Pab, 
and variation over Ab'{v)- Both formulations will be used throughout the paper. 

The idea of bipartite purifications was considered in Ref. |10i where the authors proved that every correlated 
state has, in our language, a nonzero entanglement of purification. If we would have included mixed states in the 
minimization in Eq. (|^) and used the entanglement of formation as the entanglement measure, then the defined 
quantity would be equal to the entanglement of formation of p, since the optimal extension of p is p itself. 

We put some simple bounds on Ep{p). Intuitively, 'the amount of quantum correlation in a state is smaller than or 
equal to the total amount of correlation', or Ef{p) < Ep{p). To prove this lower bound, let |V'p) = j V) A'\i) B'®\'^ij) 
be the purification that achieves the minimum in Eq. (Q). Alice and Bob locally measure the labels ia' and jb' of 
the state \tpp) such that they obtain I'f/'ij) with probability pij = (ipijlipij) . Since entanglement is nonincreasing under 
local operations, we have 

E,ip)<j:P.E(Mm)<Epip). (9) 

It is immediate that we have equality between the entanglement of formation and the entanglement of purification 
for pure states, where the optimal purification of a pure state is the pure state itself. 

An easy upper bound is Ep{p) < E{\ips){4's\) — S{pa), where pA = T^t:b{p) is the reduced density matrix in A. This 
corresponds to Ua'B' — Ia'B' or equivalently A^' = Ib' in the r.h.s. of Eq. (||) or (0). Applying the same argument 
with AA' and BB' interchanged, we obtain 

Epip)<iTnn{SipA),S{pB)), (10) 

where the purifications correspond to either completely purifying the state on A' or on B' . In general this is not the 
optimal purification, as we will see in Section 0. 

The entanglement of purification is neither convex nor concave, unlike the entanglement of formation. For instance, a 
mixture of product states, each with zero entanglement of purification, need not have zero entanglement of purification 
(for example, consider an equal mixture of |00) and |11)). On the other hand, the completely mixed state has zero 
entanglement of purification equal to zero yet it is a mixture of 4 Bell states, each with 1 ebit of entanglement of 
purification. 

Before we present continuity bounds for the entanglement of purification, we analyze the optimization problem of 
Eq. (^) in more detail. We can omit doubly stochastic maps Ab' in the optimization in Eq. (^ since they never 
decrease the entropy. Furthermore, the von Neumann entropy is concave, so that the optimum in Eq. (|^) can always 
be achieved when Ab' is an extremal TCP map. An extremal TCP map is a TCP map that cannot be expressed as 
a convex combination of other TCP maps. Choi pl has proved that an extremal TCP map with input dimension 
d has at most d operation elements in its operator-sum representation. This result will allow us to upper bound the 
dimensions of the optimal purifying Hilbert spaces, as stated in the following Lemma. 



Lemma 1 Let p act on a Hilbert space of dimension dAB = dAdB- The minimum of Eq. M) can always he achieved 

''AB 



by a state ip for which the dimension of A' is dA> = dAB OLnd the dimension of B' is dB' — d^^ (or vice versa). 



Proof. We use the formulation of the entanglement of purification as an optimization of a TCP map in Eq. (|^). 
Since the density matrix pbb' (p) is on HdB ^T^dAB : the optimal map Ab' maps TCiab into a space of some unspecified 
dimension. The optimal map Ab' can be assumed to be extremal. Theorem 5 of Choi ]TT| shows that an extremal 
TCP map A : B{Tidi) ~^ ^(^d2)[| can be written with at most di operations elements, that is, has the form 



^We have a special case when d2 = oo. The Stinespring theorem [|l^ implies that we have an operator-sum representation of 
such a map. Then Choi's results on extremality apply, bounding the number of operation elements, from which the final result 
can be proved. 
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(11) 



In our case di — (Iab- Consider implementing the TCP map by applying a unitary operation U to the input state 
with an ancilla appended. In our case, this ancilla can be taken as Alice's purifying system A' , and U acts on 
A'B' . The dimension of the ancilla A' can always be taken to be the number of operation elements. Thus we have 
dA' = dAB- The B' dimension is equal to the output dimension d2 of the optimal map A, which is unconstrained by 
the extremality condition. However, we note that the operator A(p) of Eq. has a rank of at most d^g. This is 
obtained by observing that the range of this operator is exactly that of the vectors given by all the columns of the 
matrices Vi for all i (the Vi matrices have di columns and c?2 rows). Thus, there exists a unitary operator U that 
permits the construction of a new map A' — UK whose output is confined to the first d\ dimensions of the output 
space. The operator U may be obtained explicitly via a Gram-Schmidt procedure applied to the column vectors of 
the Vi matrices. A' is also optimal, since the entropy of Eq. (^ is not changed by a unitary operation. Since the 
output space of A' has dimension d\, we conclude that ds' can be taken to he ds' = d\g. □ 

It is interesting to note that a similar minimization problem was encountered in Ref. [|l3|. There the goal was to 
use a set of noisy states for classical information transmission and we wanted to minimize the coherent information 
divided by the entropy of a quantum state under the action of a local map. 

Theorem 1 (Continuity of the Entanglement of Purification) Let p and a be two density matrices on Ti-dA ® 
Tids with Bures distance D{p,a) < e. Then 

\Ep{p) - Ep{a) I < 202?(p, a) log dAu - D{p, a) log D(p, a), (12) 

for small enough e. 

Proof: Let ji/'^) and \ipp) be the purifications of p and cr which achieve the maximum |^ in 

f^(p,a) = max |(Va|V'p)|. (13) 

Let \(j)p) and \4>a) correspond to the optimal purifications of p and cr with respect to Ep. There exists a unitary 
transformation U relating \ip'p) to \(j)p), i.e. {U (g) l)]!/']!,) — l^p)- We define the (non-optimal) purification ji/jo-) as 
{U «) l)|i/'ff) = IV'ct). Now we have 

Ep{a) - Ep{p) = i?(|0,)(0.|) - E{\cj,p){^p\) < E{\i;,){i;,\) - i?(|0,)(0,|). (14) 

We use continuity of entanglement ||l^,|l^ , Lemma (which indicates that the pure state has support on a space of 
dimension at most d^g), and the fact that |(i/'cr|0p)r— K^^lV-'p)! = F{p.a) to bound 

Ep{a) - Ep{p) < 5Dip, a) log d^ - 2D{p, a) log D{p, a). (15) 

for small enough D{p,a). We can obtain the full bound in Eq. ([l^ ) by alternatively relating \tp'^) to the optimal 
purification by a unitary transformation U . □ 

It is fairly straightforward to prove monotonicity of the entanglement of purification from monotonicity of entan- 
glement: 

Lemma 2 (Monotonicity of the Entanglement of Purification) The entanglement of purification of a density 
matrix p is nonincreasing under strictly local operations. Let Alice carry out a local TCP map Sa on the state p. We 
have 

Ep{{SA®r){p))<Ep{p). (16) 
Let Alice carry out a local measurement on p through which she obtains the state pi with probability p^. We have 

Y,P^Ep{p,)<Ep{p). (17) 

i 

Let ChOq be a local operation assisted by m qubits of communication. The entanglement of purification obeys the 
equation 



Ep{CLOq{p)) < Ep{p) + m. 



(18) 



Proof. Let \ipp) be the optimal purification of p. This optimal purification is related to some purification of 



(iSa(8i1)(/o) by a unitary transformation on Alice's system only. Then Eq. (16) follows from the fact that entanglement 



is nonincreasing under local partial traces. The state \ipi) = — A^<siB\tp) where Ai corresponds to a measurement 

^y{ip\AiA,(g,iB\->p) 

outcome of Alice, is some purification of pi. The entanglement is nonincreasing under local operations and thus 

Ep{p) = i?(|Vp)(^p|) > > Y.P^^p(P^)- (19) 

i i 

For the last inequality, let Alice and Bob start with the entangled state \ipp) and carry out their LOq protocol. By 
subadditivity of entropy, the entanglement of this state can increase by at most m bits when m qubits of communication 
are sent (back and forth). Thus the entanglement of the final state which is some purification of CLOq{p) is smaller 
than or equal to Ep{p) + m. □ 

Now we are ready to prove our main theorem: 

Theorem 2 The entanglement cost of p on TLd ® 'Ha without classical communication equals E^oqip) — E'^{p). 

Proof The inequality E^oqip) l£ ^^{p) uses entanglement dilution. Let k be the number of copies of p for which 
the regularized entanglement of purification E^ is achieved. One way of making many (p) copies of p^^ out of EPR 
pairs and o{p) < o(j)k) classical communication, is to first perform entanglement dilution on the EPR pairs so as to 
create (an approximation to) the purification and then trace over the additional registers to get p'^'^P. The other 

inequality E^{p) < E^oqip) can be proved from monotonicity and continuity of the entanglement of purification. 
We start with n EPR pairs which have Ep equal to n. The LOq process for creating an approximation pk to p'^^ 
using o(fc) qubits of communication, increases the entanglement of purification by at most o{k) bits, see Lemma H, 
or Ep{pk) < n + o(fc). Using the continuity of Theorem |l| and dividing the last inequality by k and taking the limit 
/c ^ oo gives E^{p) < ELOqip)- ° 



III. MUTUAL INFORMATION LOWER BOUNDS 



The entanglement cost E^oq is a measure of the quantum and classical correlations in a quantum state. The 
quantum and classical mutual information of a quantum state are similar measures that capture correlations in a 
quantum state. How do these measures relate to the new correlation measure? The quantum mutual information 
IqiPAs) is defined as 

Iq{pAB) = S{pA) + S{pB)-S{pAB). (20) 

We define the classical mutual information of a quantum state Ic{pab) as 

Ic{pab) = max H{pa) + H{pb) - H{pab)- (21) 

Ma-Va-Mb-Vb 

Here local measurements Ma and Mb give rise to local probability distributions pa and pb- The classical mutual 
information of a quantum state is the maximum classical mutual information that can be obtained by local measure- 
ments by Alice and Bob. Both quantum as well as classical mutual information share the important property that 
they are non-increasing under local operations (LO) by Alice and Bob. For the classical mutual information, this 
basically follows from the definition Eq. (pT|). The definition itself as a maximum over local measurements makes 
sense since the classical mutual information of a probability distribution is non-increasing under local manipulations 
of the distribution. The proof of this well known fact is analogous to the proof for the quantum mutual information 
which we will give here for completeness. 

We can write the quantum mutual information as 

IqipAB) ^ S{pab\\pa(S) pb), (22) 
where S'(.||.) is the relative entropy. The relative entropy is nonincreasing under any map A (cf. Ref. |l6[| ), i.e. 

S{k{pAB)\\^^{pA® pb)) < S{pab\\pa<»Pb)- (23) 

When A is of a local form, i.e. A^ (8) As, the l.h.s. of this equation equals the quantum mutual information of the 
state (A^ (g) Ab){pab) and thus the inequality /^((Aa Ab){pab)) < Iq{pAB) is proved. 
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A. Proof of Lower bounds 



We show that the quantities Iq{p)/2 and the regularized classical information I^{p) = lim„_^oo ^'^^ — ^ Eire both 
lower bounds for the entanglement cost E^oq- The argument is similar to the proof of the E'^ lower bound on E^oq 
in Theorem || (The reasoning is in fact a special case of Theorem 4 in Ref. (cf. Ref. applied to the class LOq 
instead of the original LOCC.) 

We start with a number, say fc, of EPR pairs which have Iq = 2k and Ic equal to k 0. In the limit of large n , 
the ratio k/n is the entanglement cost E^oqip)- We apply the LOq map C which uses o(n) communication to obtain 
an approximation p„ to p**". Since the quantum mutual information and the classical mutual information can only 
increase by o{n) by the process £ applied to the initial EPR pairs, see Lemma |3[ it follows that 

IqiPn)<o{n) + 2k, (24) 

and similarly 

IciPn) < o{n) + k. (25) 

The last step is to relate the mutual informations of pn to the mutual informations of p®". For this, we need a 
continuity result of the form 

\Iq/M)-Iq/c{p)\<Clogd\\p-a\\,+0{l). (26) 

for p,a on Tld, \\p — (t||i sufficiently small and C is some constant ^. Below we will prove these desired continuity 
results. We can divide Eqs. ( p^ ) and ( p5| ) by n and take the limit of large n. We use the continuity relation of Eq. 
( p6| ) and the fact that in the large n limit /5„ tends to p®". Thus we have 

lini = Iq{p) < 2ELOq{p), (27) 

where we used that the quantum mutual information is additive, and similarly 

I^ip) < ELOqip)- (28) 
What remains is to prove the continuity relations and the nonincrease modulo o{n) under LOq operations. 



1. Continuity of Mutual Informations 

The continuity of the quantum mutual information Iq{p) can be proved by invoking Fannes' inequality p9[ | and 
Ruskai's proof of nonincrease of the trace-distance under TCP maps [ pO[ . Let p and a be two density matrices which 
are close, i.e. ||p — ct||i = Ti\p — a\ < e for sufficiently small e. We have 

\Iq{pAB) - Iqi'TAB)] < \S{pa) ~ S{aA)\ + \S{pb) - S{aB)\ + \S{aAB) ~ S{pab)\, (29) 
which can be bounded as 

\Iq{PAB) - Iq{crAB)\ < 3\0gdAB\\PAB " CTAbHi + 5r]i\\PAB ~ CTabHi), (30) 

where r]{x) — — xloga; and — cr||i < 1/3. 

It is not hard to prove the continuity of the classical information of a quantum state, again using the nonincrease 
of ||.||i under TCP maps. Let and Mg be the optimal measurement achieving the classical mutual information 



^One can prove that Jc < fc by observing that any local measurement that is not projecting in the Schmidt basis is a noisy 
version of the measurement that does project in the Schmidt basis. In other words, the probability distribution of any set of 
local measurements can be obtained from the probability distribution of the Schmidt basis measurement by local processing, 
which does not increase the classical mutual information. 

^We can alternatively write down a continuity relation using the Bures distance. Since the trace-distance ||.||i and the Bures 
distance are equivalent distances, one continuity relation follows from the other and vice versa. 
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Ic{p)- Under this measurement the states p and a, which is, say, close to p, go to probabihty distributions p'^{i,j) 
and j) which are close again, i.e. — < ||p — cr||i. We have that 

h{a)^Icip)<I{pn-Hp'')<^ogk\y^p'^\U + 0{l), (31) 

where k is the number of joint outcomes in the optimal measurement (M'^ , Mp) and / is the classical mutual infor- 
mation of a joint probability distribution. The last inequality in Eq. ( pl| ) could in principle be derived from Fannes' 
inequality, using diagonal matrices, but it is a standard continuity result in information theory as well. To finish 
the argument, we should argue that k, the number of joint measurement outcomes is bounded. The classical mutual 
information / is a concave function of the joint probability p{i,j) Therefore only extremal measurements Ma 

and Mb need to be considered in the optimization over measurements. An extremal measurement has at most (P 
outcomes when acting on a space of dimension d [p2| and thus k < d\g. The same argument, interchanging a and p, 
can be used to upperbound Ic{p) ~ Ic{(^)- 

Lemma 3 (Monotonicity Properties of Mutual Information) Let C consist of a series of local operations as- 
sisted by m qubits of 2-way communication. The quantum mutual information obeys the inequality 

Iq{C{a))<Iq{a) + 2m, (32) 

for all states a. For the classical mutual information we have 

Ic{C{\i:){^\)<Icmm+m, (33) 

for all pure states \^p) . 

Proof. Let us first consider the quantum mutual information. We can decompose the 2-way scheme C into a 
sequence of one-way schemes. It is sufficient to prove for such a one-way scheme using m qubits of communication, 
say from Alice to Bob, that 

IgiCia)) < Igia) + 2m. (34) 

Alice's local action can consist of adding an ancilla A' in some state and apply a TCP map to the systems AA' thus 
obtaining the state ctaa'-.b- Such an action does not increase the quantum nor classical mutual information as we 
showed before. Now Alice sends system A' to Bob. We have 

IgiiAB) > I,{cjaa':b) - S{AA') + S{B) - S{AA'B) > 

S{AA') - S{A') + S{BA') - S{AA'B) > S{A) - 2S{A') + S{BA') - S{AA'B) = Iq{(TA:BA') - 2S'(A'), (35) 

where we used \S{A) — S{B)\ < S{AB) < S{A) + S{B). The quantum mutual information of the final state is 
IqicTA-.BA')- Since S{A') < to, we obtain the needed inequality. Alice could send only a part of ancilla A', but this 
does not change the bound. 

Let us now consider the classical mutual information. We may convert the entire process C into a coherent process 
C where all the measurements are deferred to the end, this does not change the amount of communication that Alice 
and Bob carry out. Thus, prior to the measurements Alice and Bob have converted the pure state into some 
pure state \(f>) whose local entropy is at most E + m where E is the entanglement of the state l^p), which is equal to 
/c(|'0)('0l) (see footnote^). Now Alice and Bob locally measure and/or trace out some registers which are operations 
that do not increase the classical mutual information. Therefore the final state £{\ip){ip\) has a classical mutual 
information that is bounded by the initial classical mutual information plus m. □. 

Remark: Note that Eq. ( |3^ ) for the quantum mutual information applies to both pure and mixed states while we 
have found mixed states that violate Eq. ( ^3|) for the classical mutual information. 

Let us state the final result once more: 

Corollary 1 ELOq{p) > Ig{p)/2 and ELOq > I^ip)- 

With this Corollary we can show that the LOq- entanglement cost of any correlated density matrix p, is nonzero ^. 
Indeed, the quantum mutual information Iq{p) of a correlated density matrix is strictly larger than zero, since S{pab) 
is strictly less than S{pa) + S{pb) (equality is only obtained when pab = Pa® Pb) and therefore ELOq{p) > 0. 

We present a simple example for which ELOqip) = ^^{p) > ^g(p)/2- 



*Note that this does not directly follow from the result in Ref. llCtl, since the entanglement of purification may be nonadditive. 
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Example 1 (All correlation is classical correlation) Consider the separable state p = J2i PiWi) i'^il ® \bi){bi\ 
where {ai\aj) = Sij and {bi\bj) = S^j . In this case Iq{p)/2 — H(p)/2. However we can show that Ep{p) > H{p). We 
have (cf. Eq. ^) p{p) = ^iPi\bi){bi\®\i){i\. Under some local TCP map K we obtain a state p! = ^^pi\bi){bi\^ Pi 
where pi are density matrices. The entropy of p' equals S{p') — '^^PiS^pi) + H(p) > H{p). The entanglement 
of purification Ep{p) may be nonadditive, so we have to consider Epi^p^"'). We have p{p'^^) = /i**" and now p' — 
j Pii • ■ ■Pin\ii ■ ■ ■ ,in){ii, ■ ■ ■ ,in\ Again the von Neumann entropy of p' is larger than or equal to 

nHija). Note that in this example we do achieve the classical mutual information lower bound. 

Here is an example where the upper and lower bounds fix the (regularized) entanglement of purification: 

Example 2 Let p be an equal mixture of the state \^o) — + |11)) "^^"^ l^i) — ^d*-"-*) ~ \^^))- Alice and Bob 

can get 1 bit of classical mutual information by both measuring in the {0, 1} basis. Thus E^oqip) ^ Ic{p) = 1, but 
ELOqip) < S{pa) < I, Eq. Therefore E^Oq = 1- 



IV. OTHER CORRELATION MEASURES: THE LOCALLY INDUCED HOLEVO INFORMATION 

In Ref. [Q the authors considered the locally induced Holevo information as a measure of classical correlations in 
the state. It is defined either with respect to Alice's measurement {Ca) or Bob's measurement {Cb) 

Ca/b{p) - iTiax SiY.pf'^'pf"') - J2pf^^S{pf^\ (36) 

Ma /Mb ^—^ 

where Ma {Mb) on p gives reduced density matrices pf {pf) with probabihty pf [pf). The classical mutual 
information I^{p) will in general be less than these quantities, since to achieve the Holevo information one may have 
to do coding. In Ref. it was shown that Ca/b are nonincreasing under local operations. We leave it as an exercise 
for the reader to prove continuity and nonincrease modulo o(n) under LOq operations (applied to some pure state), 
thus showing that the regularized versions of these two quantities are also lower bounds for E^oq- 



A. Bell-diagonal states 

We show that for Bell-diagonal states pBeii the quantity Ca (equal to Cb by symmetry of the Bell-diagonal states) 
is equal to the classical capacity of the corresponding qubit channels. By the previous arguments this give us some 
lower bounds on the regularized entanglement of purification of these states. The Bell-diagonal states are of the 
following form 

i 

where ^'o...3 are the four Bell states where |^o) is :^(|00) + |11)). The corresponding channel, -the so called generalized 
depolarizing channel-, or Pauli channel, is of the form 

Ap(-) = (38) 

i 

where ctq = 1, and o'i_2,3 are the three Pauli matrices. It is known |Q that all two qubit states with maximally 
mixed subsystems are Bell-diagonal, up to a unitary transformation Ua <8) Ub- From the isomorphism between states 
and channels ||2^,|ll|,^, it follows that all unital channels are of the form (|3|) [cf. [^], up to unitary transformations 
applied before and after the action of the channel. The classical 1-shot capacity of the quantum channel A is given 
by [160 

Ci(A) = sup x({ q^,Mp^)}), (39) 

{q^,P^} 

where x is the Holevo function of the ensemble 

X({ qi,Pi}) = sC^qiPi) -^qiS{pi). (40) 

i i 
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The optimal states pi that achieve the capacity Ci are always pure states, moreover it can be shown that the 
ensemble {qi, that achieves Ci for unital 1-qubit channels satisfies 



Y.q^\4'^){^^\^\l. (41) 



2 

Let us argue that Ca{p) = C'i(A) for a Bell-diagonal state pBeii = (lA'8)Ap)(|^o)(^o|)- Ahce's POVM measurement 
on this state commutes with the channel A^. By doing a measurement on |^'o) she can create any pure-state-ensemble 
on system i3, obeying the relation Eq. (41). This ensemble is then sent through the channel Ap. If the ensemble is 
optimal for Ci , then its Holevo information Yequals Ci and thus Ca = Ci . 

For unital 1-qubit channels Ci is given by p5| , p8[ 

C7i(A) = l-min5(A(|V')(V'|))- (42) 

We can perform the minization in the last inequality and we obtain the following formula for the capacity of a Pauli 
channel or the induced Holevo information of the Bell-diagonal states 

CA{pBeii)^Ci{Ap) = l-H{l-X), (43) 

where A is the sum of the two largest probabilities pi and H{.) is the binary entropy function H{x) = — xloga; — (1 — 
x) log(l — x). For a two-qubit Werner states of the form 

3 

pw = e|vI'o)(vI'o| + (1 - e)/3^ |vl',;)(*,|, (44) 

i=l 

we obtain 

CA = l-i?(^) foree[il], 

Ca = 1-h(^^^^ foree[0,i]. (45) 

It was shown by King that the classical capacity of unital 1-qubit channels is equal to the one shot capacity, 
or Ci = = hm„^oo iCi(p»"). Therefore Ca = = Ci, which is a lower bound on ELOq- 



V. WERNER STATES 

A numerical minimization based on Eq. (|^) was performed for the Werner states Eq. (^^ for Ep. We plot the 
results as a function of the |4'o) eigenvalue e in Fig. |l|. We permitted various output dimensions; The two curves shown 
have dim(A') — dim(i3') = 2 and dim(A') ~ dim(i?') = 4. In the first case, the initial variable of the minimization 
was determined by a random 4x4 unitary Ua'B' picked according to the Haar measure. In the second case, the initial 
point was determined by a random 16 x 4 isometry picked according to a parameterization derived from Ref. [ p9[ . 
We did not explore the largest dimensions permitted by Lemma |l], which would have required an optimization over 
a 64 X 4 isometry. 
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Eigenvalue of \'io) 



FIG. 1. Numerical bounds on Ep for Werner states. In the upper curve we restrict to dim(j4') = dim(_B') — 2; for the next 
curve, we permit dim{A') — dim(_B') — 4. The inset shows the curious beha vior o f Ep around the point where the eigenvalue 



of jv^o) approaches zero. The dotted curve is the Ca lower bound of Sec. IVA. The dashed curve is the entanglement of 



formation lower bound which vanishes when the eigenvalue is smaller than or equal to 1/2. 



It is evident from the numerics presented in the figure that the Ca bound of Eq. (^) is not achieved for the Werner 
states: the Ca lower bound is only tight at the trivial points e = 1/4 and e — 1. Our results indicate that Ep is a very 
complex function, neither concave nor convex, with several distinct regimes. In fact, we find four different regimes in 
our numerics: I) In this regime the standard purification of Eq. (0) appears to be optimal, so the U of Eq. (^) is the 
identity, and the purifying dimensions are dim(A') — 1 and dim(_B^) — 4. This regime only extends over a tiny range, 
approximately < e < 0.005. II) In the range 0.005 < e < 0.25 we find an optimal purification of the form 




In this region the Ep curve is given by Ep — — xlogx — (1 — x) log((l — a;)/3), with cc = (1 + 2e — 2-\/3\/e(l — e))/12. 
Here the purifying dimensions are dim(A') = 2 and dim(i?') = 2. Of course Ep drops to zero for the completely mixed 
state at e = 1/4. Ill) In the range 0.25 < e < 0.69 we also find purifying dimensions dim(^') = 2 and dim(i3') — 2, 
but we were unable to determine the analytical form of the purifying state or of Ep. IV) In the range 0.69 < e < 1 
the purifying dimensions were dim(A') = 2 and dim(_B') = 3. Again, we were unable to come to any analytical 
understanding of the result. Of course, Ep = 1 for e = 1, corresponding to the pure maximally entangled state. 



VI. CONCLUSION 



We have shown that the entanglement cost ELOqip) is equal to the regularized entanglement of purification. It is 
an open question whether the entanglement of purification is additive: 

Ep{p(it>p) = Ep{p) + Ep{p). (47) 
In the alternative formulation using the state ^{p) the additivity question is the following. Is the minimum in 

Tn\nS{{lAB ® ^CD){^J.AC ® ^J■BD), (48) 

AcD 

achieved by a TCP map Ac d = S ®S1 This problem is similar again to the additivity question encountered in Ref. 

where a local map could possibly lower the ratio of the coherent information and the entropy of many copies of 
a state together. 

It is interesting not only to ask the formation question with respect to this class LOq, but also consider 'the 
distillation' question. One can consider different versions. For example, how much entanglement can we distill from 
p using o(n) communication? One would expect that this quantity DLOq{p) is always zero for states for which the 



10 



entanglement cost Ec (using LOCC) is lower than the distillable entanglement D. We do not have a proof of this 
statement, relating irreversibility to a need for classical communication. 

Instead of trying to convert the correlations in p back to entanglement, we may ask what classical correlations Alice 
and Bob can establish using p. We could allow Alice and Bob to perform an asymptotically vanishing amount of 
communication in this extraction process. A little bit of communication could potentially increase the classical mutual 
information in a quantum state by a large amount (when the classical correlation is initially 'hidden'), thus this may 
not be the best problem to pose. Researchers |Q |3l] have investigated the possibly more interesting problem of 
the secret key K that Alice and Bob can establish given p where one allows arbitrary public classical communication 
between the parties. There is again more than one version of this problem, one in which Eve possesses the purification 
of the density matrix ]30|] and a situation in which Eve is initially uncorrelated with the density matrix. In Ref. ]3^ ] 
a general framework is developed to address these issues also in the multipartite setting. 

Quite recently, entanglement properties of bipartite density matrices were studied by looking at mixed extensions of 
the density matrix . It would be interesting to explore the connection between our results here on the entanglement 
of purification and this other approach. 
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